An adaptive routing strategy for packet delivery in complex networks 
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We present an efficient routing approach for delivering packets in complex networks. On delivering a message 
from a node to a destination, a node forwards the message to a neighbor by estimating the waiting time along 
the shortest path from each of its neighbors to the destination. This projected waiting time is dynamical in 
nature and the path through which a message is delivered would be adapted to the distribution of messages in 
the network. Implementing the approach on scale-free networks, we show that the present approach performs 
better than the shortest-path approach and another approach that takes into account of the waiting time only at 
the neighboring nodes. Key features in numerical results are explained by a mean field theory. The approach has 
the merit that messages are distributed among the nodes according to the capabilities of the nodes in handling 
messages. 

PACS numbers: 89.75.Fb,89.20.-a,05.70.Jk 
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I. INTRODUCTION 

The problem of traffic congestions in communication net- 
works is undoubtedly an important issue. The problem is re- 
lated to the geometry of the underlying network, the rate that 
messages are generated and delivered, and the routing strat- 
egy. Many studies have been focused on spatial structures 
such as regular lattices and the Cayley tree QjJ, |2j, Ij, U, |5l |6[ |7j, 
• Random networks and scale-free (SF) networks have 
also been widely studied. The former is homogeneous with a 
Poisson degree distribution; while the latter typically exhibits 
a power-law degree distribution of the form P(k) ~ k^ 1 sig- 
nifying the existence of nodes with large degrees. SF net- 
works are found in many real-world networks, such as the 
Internet World Wide Web (WWW), and metabolic network 
11311 . A standard model of SF networks is the Barabasi 
and Albert (BA) model of growing networks with preferential 
attachments lfl4l [l5tl . The BA model gives a de gree distribu- 
tion of P(k) ~ k~ 3 and is non-assortative 1161. 11711 . i.e., the 
chance of two nodes being connected is independent of the de- 
grees of the nodes concerned. While there are other variations 
on the BA models that giv e a degree exponent that deviates 

from 3 EE 113 EH III, 

the BA model still serves as the 
basic model for SF networks. In the present work, we study 
how a dynamical and adaptive routing strategy would enhance 
the performance in delivering messages over a BA scale-free 
network. 

In communication models, the nodes are taken to be both 
hosts and routers, and the links serve as possible path- 
ways through which messages or packets are forwarded to 
their destination. Early studies assumed a constant (degree- 
independent) packet generation rate A and a constant rate of 
delivering one packet per time step at each node. As A in- 
creases, the traffic goes from a free-flow phase to a congested 
or jamming phase. Obviously, such models are too simple for 
real-world networks. More realistic models should incorpo- 
rate the fact that nodes with higher degrees would have higher 
capability of handling information packets and at the same 



time generate more packets. Models with degree-dependent 
packet-delivery rate of the form (l+0kj) B23ll24ll and degree- 
dependent packet-generating rate Il24ll of the form Afc; have re- 
cently been proposed and studied, with routing strategy based 
on forwarding messages through the shortest path to their des- 
tination. Implementing this routing strategy in SF and random 
networks indicate that it is easier to lead to congestions in SF 
networks than in random networks J23l l24ll . It is because the 
nodes with large degrees in SF networks are on many shortest- 
paths between two arbitrarily chosen nodes, i.e., large be- 
tweenness [25]. Many packets will be passing by and queue- 
ing up at these nodes enroute to their destination. In a random 
network, jamming is harder to occur as the packets tend to be 
distributed quite uniformly to each node. 

A good routing algorithm is essential for sustaining the 
proper functioning of a network 1261, l27ll28l, 12911 . The shortest- 
path routing approach is based on static information, i.e., once 
the network is constructed, the shortest-paths are fixed. To im- 
prove routing efficiency, Echenique et al. IBOlblll proposed an 
approach in which a node would choose a neighboring node 
to deliver a packet by considering the shortest-path from the 
neighboring node to the destination and the waiting time at the 
neighboring node. The waiting time depends on the number 
of packets in the queue at a neighboring node at the time of de- 
cision and thus corresponds to a dynamical or time-dependent 
information. This algorithm performs better than the shortest- 
path approach, as packets may be delivered not necessarily 
through the shortest-path and thus the loading at the higher 
degree nodes in a SF network is reduced. The approach has 
also been applied to networks with degree-dependent packet 
generation rate [32]. Recently, Wang et al. lO proposed an 
algorithm that tends to spread the packets evenly to nodes by 
considering information on nearest neighbors. However, the 
delivering time turns out to be much longer than that in the 
shortest-path approach as the packets tend to wander around 
the network. 

In the present work, we propose an efficient routing strategy 
that is based on the projected waiting time along the shortest- 
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path from a neighboring node to the destination. The algo- 
rithm is implemented in BA scale-free networks, with degree- 
dependent packet generating and delivering rates. Results 
show that jamming is harder to occur using the present strat- 
egy, when compared with both the shortest-path approach and 
the Echenique's approach. Key features observed in numeri- 
cal results are explained within a mean field treatment. The 
present approach has the advantage of spreading the packets 
among the nodes according to the degrees of the nodes. In this 
way, every node can contribute to the packet delivery process. 

The paper is organized as follows. The model, including 
the underlying network, the packet generation and delivery 
mechanisms, and routing strategy, is introduced in Sec. II. In 
Sec. Ill, we present numerical results and compared them with 
those of the other routing strategies. We also explain key fea- 
tures within a mean field theory. We summarize our results in 
Sec.IV. 



II. MODEL 

The underlying network structure is taken to be the 
Barabasi-Albert (BA) scale-free growing network with N 
nodes 11511 . Starting with too nodes, each new node enter- 
ing the network is allowed to establish to new links to exist- 
ing nodes. Preferential attachment whereby an existing node 
i with a higher degree ki has a higher probability II; ~ fej to 
attract a new link is imposed. The mean degree of the net- 
work is (k) = 2m and the degree distribution P(k) follows a 
power-law behavior of the form P(k) ~ fc~ 3 . 

The dynamics of packet generation and delivery is imple- 
mented as follows. Due to the inhomogeneous nature of the 
BA network, it is more natural to impose a packet genera- 
tion rate that is proportional to the degree of a node. At each 
time step, a node i creates Xki new packets. The fractional 
part of Xki is implemented probabilistically. A destination is 
randomly assigned to each created packet. The newly cre- 
ated packets will be put in a queue at the node and delivery 
will be made on the first-in-first-out basis. The packets in the 
queue may consist of those which are created at previous time 
steps and received from neighboring nodes enroute to their 
destination. We also assume a packet delivery rate that is pro- 
portional to the degree of a node B2411 . At each time step, a 
node i delivers at most (1 + (3ki) packets to its neighbors. The 
fractional part of (1 + (3ki) is implemented probabilistically. 
A larger (3 implies a higher packet-handling capability, but it 
would translate into higher cost or capital. Here, the param- 
eters A and (3 are taken to be node-independent. A packet is 
removed from the system upon arrival at its destination. For 
a given generation rate characterized by A, there exists a crit- 
ical value of the delivery rate j3 c such that for delivery rates 
(3 < (3 C , packets tend to accumulate in the network resulting 
in a jamming phase; while for (3 > (3 C , a non-jamming phase 
results as there are as many packets delivered to their destina- 
tion as created. A better performance is thus characterized by 
a smaller value of j3 c . 

The novel feature of the present work is the routing strategy 
or the selection of a neighbor in delivering a packet. The idea 



is to choose a neighbor that would give the shortest time, in- 
cluding waiting time, to deliver the packets along the shortest 
path from the chosen neighbor to the destination. Consider a 
packet with destination node j leaving node i. Each of the fc; 
neighbors of node i has a shortest path to the destination node 
j. The shortest path refers to the smallest number of links 
from a node to another. However, due to the possible accu- 
mulation of packets at each node, the number of time steps it 
takes to deliver the message may be different from the number 
of links along the shortest path. Consider a neighbor labelled 
£ of the node i. We label the shortest path from node I to j 
by {SP : £,j}. Along this path, we evaluate the following 
quantity for the node £: 

- :-/':••...: ' 

where the sum is over the nodes along the shortest path 
{SP : £, j}, excluding the destination. Here, n s is the number 
of packets accumulated at node s, at the moment of decision. 
Thus, d(£) is an estimate of the time that a packet would take 
to go from node I to the destination j through the shortest 
path. Node i would choose a neighboring node with the mini- 
mum d(£) to forward the packet, i.e., the selection is based on 
min{d(£), £ £ {i}}, where {i} is the set of ki nodes consist- 
ing of the neighbors of node i. This procedure is repeated for 
each node and each packet in every time step. For a network 
far from jamming, each node can handle all the packets in ev- 
ery time step. In this free-flow situation, the quantity d(£) sim- 
ply measures the shortest path dg.j from £ to j. When packets 
are queueing up at the nodes, however, a delivery mechanism 
based on d(£) takes into account of the queueing time and may 
not pass the packet to a neighboring node that is closest to the 
destination. 

To justify our routing scheme, we will compare results with 
two other routing strategies widely studied in the literature. 
Using the same packet generating mechanism, the shortest- 
path approach selects a neighbor with the shortest path to the 
destination for forwarding a packet. Echenique et al. ll30l[3Tll 
proposed an approach that takes into account of the waiting 
time. For a delivering rate of one packet per time step, they 
proposed to choose a neighbor that has a minimum value of 
hdej + (1 — h)ri£, where dg j is the shortest path length from 
node £ to j. The parameter h is a weighing factor, which can 
be taken as a variational parameter and h w 0.8 is found to 
give the best performance. The Echenique's approach thus 
accounts for the waiting time only at the neighboring nodes. 
For a delivery rate of (1 + (3ki), a modified Echenique's ap- 
proach is to choose a neighboring node with a minimum value 
of 

5 e = hd itj + (l-h)-^—. (2) 

We have checked that for a given value of A, the smallest value 
of f3 c is attained for values of h ~ 0.8 to 0.85. In what follows, 
we will use a value of h — 0.8 for the Echenique's approach 
given by Eq.(2). 
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FIG. 2: Numerical results for the quantity r\ as function of /3 for 
FIG. 1: (a) Numerical results for the average number of packets per A = 0.02. The critical value /3 C separates the behavior of n = for 
node < n(t) > and (b) the average delivering time < T > as a P > (3 C and j] for (3 < j3 c . 
function of time for A = 0.02. Lines from top to bottom correspond 
to (3 = 0.04, 0.048 and 0.07. 
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III. RESULTS AND DISCUSSION 

The different phases in a network can be illustrated by look- 
ing at the average number of packets per node at a given time 
(n(t)) and the average time for a packet to remain in the net- 
work or the delivering time (T). We take mo = 3 and m — 3 
and construct a BA scale-free network of N = 1000 nodes. 
Figure Q] shows the results of (n(t)) and (T) as a function of 
time for a fixed value of A = 0.02. As (3 increases, there are 
distinct behavior. For values of (3 smaller than some critical 
value (3 C (X), (n(t)) grows almost linearly with time after the 
transient (see FigQ~|a)). This corresponds to a jamming phase. 
As j3 increases, the slope in the long time behavior decreases, 
indicating a slower accumulation of packets in the network as 
the ability of handling packets (3 increases. For (3 > /3 C (A), 
(n(i)) becomes independent of time in the long time limit. 
This corresponds to a non-jamming phase. Similarly behav- 
ior is exhibited in (T). In the jamming phase, (T) increases 
with time monotonically, due to the increasing waiting time in 
the queues at intermediate nodes as a packet is forwarded to 
its destination. Fewer packets are delivered to their destina- 
tion than generated. In the non-jamming phase, (T) becomes 
independent of time in the long time limit. In this regime, fur- 
ther increasing (3 will lead to smaller (n(t)) and shorter (T) 
in the long time limit until these quantities saturate. This is 
possible since a non-jamming phase corresponds either to the 
case in which all the packets at the nodes are forwarded ev- 
ery time step or steady queues of packets exist at the nodes. 
In both cases, the number of packets does not increase in the 
long time limit. The former case is the free-flow phase, while 
the latter is reminiscent of the synchronized phase in vehicular 
traffic flows [34] in which the packets undergo a stop-and-go 
behavior. For (3 = 0.07 > (3 C for example, (T) ks 9.5, which 
is somewhat larger than the average shortest distance or diam- 
eter D w 3.332 of the network. This indicates that, due to 
the routing strategy in forwarding a packet, the dynamics in 
the free-flow phase is different from that of the shortest-path 
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FIG. 3: The critical value /3 C (A) for three different routing ap- 
proaches for forwarding packets: the present approach (circles), the 
Echenique's approach with h = 0.8 (stars), and the shortest-path 
approach (squares). The lines are guides to eye. 



approach. 

The critical value /3 C (A) can be determined by considering 
the quantity 



r) = lim 



1 < An > 



t->°o 2mA At 



(3) 



where An = n(t + At) — n(t) and the average is over all 
the nodes at a time t. This quantity rj G [0, 1] is basically the 
slope of (n(t)) in the long time limit. In the non-jamming 
phase, the slope vanishes and n = 0; while in the jamming 
phase, n > 0. Figure |2] shows 77 as a function of (3, for a fixed 
value of A = 0.02. The critical value (3 C can be identified as 
the value that separates the 77 = and 77 ^ behavior. We 
carried out similar calculations for different values of A and 
determined /3 C (A). The results are shown in FigO (circles). 
We will explain the form of (3 C (\) using a mean field theory. 
The curve (3 C (X) can also be regarded as a phase boundary in 
the \-(3 plane, separating the jamming phase below the curve 
and the non-jamming phase above the curve. 

To show the superior performance of our routing strategy, 
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FIG. 4: The average number of packets n k at a node of degree k 
as a function of fc for A = 0.02 for (a) (3 = 0.06(> f3 c ) and (b) 
/3 = 0.04(< /3 C ). In each case, results are shown at three different 
times of t = 100 (circles), 200 (stars) and 300 (squares) time steps. 



we also performed calculations using the shortest-path ap- 
proach and the Echenique's approach with h = 0.8 in Eq.(2). 
The same degree-dependent packet generating mechanism is 
used. Results of (3 C (X) for these two models are also shown 
in Fig|3] for comparison. The present approach gives the 
best performance. For a given A, we see the improvement 
in performance from the shortest-path approach through the 
Echenique's approach to the present approach, signified by 
the drop of f3 c . For the shortest-path approach, it has been 
shown Il24ll that /3 C (A) follows the functional form of 



I3 S C P = aD(X - A 



sp ) 

in ! n ) 



(4) 



where a 2 and A^„ — l/(aDk max ) with D being the di- 
ameter and k max the maximum degree of the network. For 
A < A^„, /3f p = 0. With the present approach, C (X) 
follows a similar functional form, but with a higher value of 
\ m in and a smaller prefactor that gives the slope. Both the 
present approach and the Echenique's approach perform bet- 
ter than the shortest-approach approach because packets are 
re-directed to other nodes when there are long queues at the 
hubs. 

The better performance of the present approach is achieved 
by spreading the packets among the nodes so that the number 
of packets at a node is proportional to the degree fc of the 
node in the free-flow phase. We use a mean field approach to 
illustrate this point. Let n k be the average number of packets 
at the nodes with degree fc. In the free-flow phase where n k < 
1 + j3k, we have 



dn k (t) 
dt 



= Xk—n k (t)+k 



k' = k mi 



P(k'\k)^P--X(k). (5) 



The first and second terms denote the packets generated at the 
node and delivered to neighboring nodes, respectively. The 
third term accounts for the packets delivered into the node 
from its neighboring nodes. Here P(fc'|fc) is the conditional 
probability that a node of degree fc has a neighbor of degree 
fc' and the sum runs from the minimum degree fc„ 



in the network. In the free-flow regime, the packets that are 
removed upon arrival at their destination can be assumed to be 
fc-independent and approximated by the term A(fc). The non- 
assortative feature of BA networks llfj, Il7ll gives P(fc'|fc) = 
k'P(k') I (fc), where P(fc') is the degree distribution. After the 
transient behavior, dn k /dt = and we have 



^ < n > \ , w> X 

»fc = A+— — fc-A fc , 
< fc > 



(6) 



where (n) = Y/k^k - p ( k ') n k> is the mean number of 
packets per node. Thus for fc > (fc) = 6, n k ~ fc in the 
free-flow phase after the transient. Figure HJa) shows the nu- 
merical results obtained by averaging the number of packets 
on the nodes with degree fc at different times (time t = 100, 
200, 300 time steps) of a run. In the free-flow phase, n k ~ fc 
and becomes time-independent after the transient, as shown in 
Figj4ja) for the case of j3 = 0.06 and A = 0.02. This behavior 
is consistent with that in Eq.©. 

For the jamming phase, numerical results (see FigHfb)) 
show that (i) n k ~ fc at a fixed instant and (ii) n k increases 
with time for fixed value of fc. This behavior can be under- 
stood provided that the packets are still distributed among the 
nodes in proportion to the degree fc of a node via our strategy. 
In this phase, the long time behavior is characterized by an 
increasing accumulation of packets and the delivery to desti- 
nations becomes negligible compared with packet generation. 
With rife > 1 + (3k for all nodes and ignoring the removal of 
packets, Eq.(6) is modified to 



dn k (t) 
dt 



= Afc-(l + /?fc) + fc 2 P(k'\k) 1 + f ^ 



k '=k min 



= -l + fc(A + 



< k >' 



It follows that n k (t) increases with time t as 



n k (t) =n fe (0) + (fc(A + 



< fc > 



) - 1)*, 



(7) 



(8) 



which describes very well the features in Fig|4fb). Thus, the 
present approach has the effect of reducing (increasing) the 
probability of passing packets to neighbors with high (low) 
degrees when there are long (no or short) queues, resulting in 
a distribution of packets according to the degrees of the nodes. 

A rough estimate of f3 c (A) can be obtained by equating 
n k in the free-flow phase to (1 + (3 c k). In particular, taking 
n k ma * = 1 + Pckmax, we get from Eq.© that 



lie = A- A 



A - A 



(fc) < n > 1 



(fc) , (fl-l)EAi, 1 



D- 



(k) 



X- 



(D - (fc)/fc 

'max J^n 



d x- 



Dkr, 



(9) 
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where D is the average number of nodes that a packet passes 
through from its origin to the destination, which is the diame- 
ter of the network in the free-flow phase. The last line is valid 
for k max 3> (k). Comparing with Eq.© for the shortest-path 
approach, we note that X min = 1/ (Dk max ) > A^f n and the 
prefactor D, which gives the slope in Fig[3] is smaller than 
that in the shortest-path approach. These features are con- 
sistent with numerical results. In particular, for N — 1000 
nodes, we found that D w 3.332 and k max « 85, giving 
A m m ~ 0.007, which is in reasonable agreement with numer- 
ical results in Figf3] 

IV. SUMMARY 

In summary, we have proposed an efficient routing strategy 
on forwarding packets in a scale-free network. The strategy 
accounts not only for the physical separation from the des- 
tination but also on the waiting time along possible paths. 
We showed that our strategy performs better than both the 
shortest-path approach and the Echenique's approach. Ana- 
lytically, we construct a mean field treatment which gives re- 
sults in agreement with observed features in numerical results. 
Our routing strategy has the merit of distributing the packets 
among the nodes according to the degree, and hence handling 
capability, of the nodes. Although our discussion was carried 
out on BA networks, we believe that our approach is also ap- 
plicable in other spatial structures. 

We end by comparing the three different routing strategies 
in more general terms. The shortest-path approach depends 
entirely on geometrical information that is static. Once the 



origin and the destination of a packet is known, the shortest- 
path is fixed. This strategy is non-adaptive, i.e., it will not be 
change with time. The Echenique's approach considers both 
geometrical and local dynamical information. By considering 
the waiting time at a neighboring node, a packet from a node i 
to a destination j will not always follow the same path. Thus, 
the Echenique's approach is a strategy that is adaptive, i.e., a 
decision based on the current situation. The present strategy, 
like the Echenique's approach, is also adaptive and makes use 
of global information in which all the waiting times along a 
path are taken into consideration. We see that by allowing 
for adaptive strategies and taking more information into con- 
sideration, a better performance results. This line of thought 
is in accordance with that in complex adaptive systems D5I1 
whereby active agents may adapt, interact, and learn from past 
experience. It should be, however, noted that it pays to be bet- 
ter. The shortest-path approach does not require update of the 
routing strategy. The Echenique's approach and the present 
approach require continuing update of the number of packets 
accumulated at the nodes. Such updating plays the role of a 
cost, with the payoff being the better performance. Practical 
implementation would have to consider the balance between 
the cost and the payoff. 
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